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It is well understood that the studies of correlations between produced 
particles, the effects of coherence and chaoticity, an estimation of parti- 
cle emitting source size play an important role in high energy physics [1]. 
First of all, we mean the investigation of the space-time extension or even 
squeezing of particle sources via the multiparticle quantum-statistics cor- 
relation. We consider the two-photon correlation function that can pro- 
vide the space-time information about the Higgs-boson source in thermal 
environment and estimate the Higgs-boson mass for the first time. 



1. Very soon the Large Hadron Collider (LHC) at CERN will provide 
particle physicists with treasure of data. These data allowed a precise 
measurement of many important parameters in order to test their con- 
sistency and to hunt the Higgs boson. At the same time, little is known 
about Higgs-gauge bosons interplay in particular in theoretical aspects. 
The effect of Bose-Einstein correlations (BEC) is clear and undeniable 
part of this theory, complicating the quantum statistical description of 
multi-lepton final states. It is known that the wave function describing a 
system of two identical bosons should be symmetric under permutation 
of these particles. As a consequence, the four-momentum differences be- 
tween two bosons will be smaller than in the world where Bose-Einstein 
statistics would not apply. 

In Higgs searching program at the LHC we suppose that two photons 
produce mainly through Higgs-boson decay. ^From the theoretical point 
of view the decay of Higgs into two photons emerges on the one-loop 
level: Higgs-boson transferred in pair of a quark and antiquark or vec- 
tor bosons, or even in scalar particles which annihilate in two photons. 
Experimentally the 7r° background constitutes a major difficulty which, 
however, to some extent, can be taken care of by measuring the photon 
pair invariant mass. 



The importance of two-photon correlations investigation in Higgs physics 
relates to the fact that the Higgs-boson mass is estimated via the correla- 
tion radius which defines the geometrical size of the two-photon source. 
Hence, the two-photon correlation function can be strongly dependent 
on the space-time properties of the Higgs decays. 

In the quark-loop scheme the non-relativistic bremsstrahlung formula 
for the current in quark-antiquark interaction for one photon emission 
(with four-momentum k) is 

J X (k) = ^p-e x (k)e- k °^, (1) 

where p is the difference between the momenta of a quark and antiquark 
with the mass m q , £\{k) is the vector of linear polarization of a photon; 
e means the phenomenological parameter which depends on the initial 
energy. For two sources the transition current is 

2 

j\k) = J2 eikVn in(k). (2) 

n=l 

Index n labels the independent quark-antiquark interactions taking place 
at different space-time points y n which are considered to be randomly dis- 
tributed in the space-time volume of the (photon) source. We suppose 
the gamma-quanta are created by some random currents (sources) J^(x) 
through the Lagrangian density L int = —gJ fJi (x)A^(x). The neutral 
weak current of charged leptons is J l Jx) = —g l(x) 7 M l(x) with I : e, /i, r, 
for which the problem can be solved exactly. The currents J^x) and 
Jpix) exist in a restricted space-time region and to be chaotic and ran- 
domly disturbed by external fields (forces). The restricted domain is 
characterized by the internal stochastic scale L st the meaning of which 
will be explained later. 

By studying BEC of identical particles (e.g., photons ), it is possible 
to determine the time scale and spatial region over which particles do not 
have the interactions. Such a surface is called as decoupling one. In fact, 
for an evolving system such as pp collisions, it is not really a surface, 
since at each time there is a spread out surface due to fluctuations in 
the last interactions, and the shape of this surface evolve even in time. 
The particle source is not approximately constant because of energy- 
momentum conservation constraint. 
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Actually, the second order distribution function N 12 (ki, k 2 ) = (\J Xl (k 1 ) J X2 (k 2 )\ 2 } 
normalized to Ni(ki) ■ N 2 {k 2 ) with one-particle distribution functions 
Niiki) = (\J Xl (h)\ 2 } (i = 1,2) formally defines the probability to find 
two photons with momenta k\ and k 2 issued at y\ and y 2 . The crossing 
momenta has to be taken into account. 

In this work, we make an attempt to demonstrate that the problem 
of properties of the genuine interactions can be explored using experi- 
mental data which can be collected by CMS detector at the LHC. These 
data can be analyzed through the compared measures of some inclusive 
distributions and final state correlations. 

One of the aims of this paper is to carry out the proposal for the 
experimental measurements of virtual 77 pair correlations. 

This exploration is theoretically supported by the quantum field the- 
ory model approaches [2-9] at finite temperature, (QFTp), where one of 
the main parameters is the temperature of the particle (emitting) source 
under the random force operator influence. 

We propose that photons do not strongly interact with (produced) 
medium: they carry information about early stage of reaction. Our prag- 
matic definition is: photons are produced not from hadronic decays. Any 
source of real gamma-quanta produces virtual photons with very low 
mass. If the momentum of the virtual gamma-quantum is sufficiently 
small, the source strength should be small as well. The real gamma- 
quantum can be measured from the virtual yield which is observed as 
low mass of the spectrum for lepton-antilepton pair. 

The main channels are the two-photon production pp — > Higgs — > 
7*7* — > 2/i~2/i + , 2e~2e + , e~e + /i~/i + , ... in pp collisions. An efficient 
selection of leptons needs to be performed according to the following 
criteria. First, all leptons were required to lie in the pseudorapidity 
range covered by, e.g., the CMS muon system that is, < 2.4. Second, 
the leptons were required to be unlikely charged in pairs. 

The di-lepton channel is especially promising from the experimental 
point of view, since it is expected that the experimental facilities related 
for LHC will make it possible to record muons of energy in the TeV 
range with a resolution of about a few percent and an efficiency close to 
100 %. Moreover, this channel is characterized by a maximum signal-to- 
background ratio in the energy region being considered. 

2. A pair of identical bosons with momenta p\ and p 2 and the mass 
m produced incoherently from an extended source will have an enhanced 
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probability C 2 (pi,j? 2 ) = N 12 (p 1 ,p 2 )/[N 1 (p 1 ) ■ N 2 (p 2 )] to be measured in 
terms of differential cross section er, where 

to be found close in 4-momentum space when detected simultaneously, 
as compared to if they are detected separately with 

1 da d s pi 



NM = — — , dVLi = - — -p , E v . = Jpf + m 2 , i = 1, 2. (4) 

In an experiment, one can account the inclusive density p2{Pi,P2) which 
describes the distribution of two particles in Q (the sub-volume of the 
phase space) irrespective of the presence of any other particles 

, ,11 d 2 n 2 

P2{P1,P2> 



2! n events dpi dp 2 ' 

where n 2 is the number of particles counted in a phase space domain 
(pi + dpi,p 2 + dp 2 ). The multiplicity N normalizations stand as 

p(p)dp=(N), 



PziPu&i) dpi dp 2 = (N(N - 5 12 )), 



n 



where (N) is the averaged number of produced particles. Here, 612 = 
for different particles, while 5u = 1 in case of identical ones (coming from 
the same event). 

On the other hand, the following relation can be used to retrieve the 
BEC function C 2 (Q): 

<*(«) - <»> 

where N(Q) in general case refer to the numbers of off-shell photons with 
Q = \J-(pi-P2)^(pi-P2Y = VM 2 -4m 2 . (6) 

In definitions (jSJ) and ([H]), N re f is the number of particle pairs without 
BEC and p^ = (u)i,Pi) are four-momenta of produced photons (i = 1, 2); 

M = \ {Pi + P2H is the invariant mass of the pair of photons. 



An essential problem in extracting the correlation is the estimation 
of the reference distribution N re f(Q) in Eq. (J5J). If there are other corre- 
lations beside the Bose-Einstein effect, the distribution N re ^ (Q) should 
be replaced by a reference distribution corresponding to the two-particle 
distribution in a geometry without BEC. Hence, the expression (j3]) rep- 
resents the ratio between the number of 77 pairs N(Q) in the real world 
and the reference sample N re f(Q) in the imaginary world. Note that the 
reference sample can not be directly observed in an experiment. Differ- 
ent methods are usually applied for the construction of reference samples 
[1], however all of them have strong restrictions. One of the preferable 
methods is to construct the reference samples directly from data. For 
our aim for reference sample N re f(Q), it is suitable to use the pairs 77 
from different (mixed) events. 

It is commonly assumed that the maximum of two-particle BEC func- 
tion C2{Q) is 2 for pi = P2 if no any distortion and final state interactions 
are taking into account. 

In general, the shape of C 2 (Q) is model dependent. The most simple 
form of Goldhaber-like parameterization for C 2 (Q) [10,11] can be used 
for experimental data fitting in particle physics (see, e.g., [1] and the 
references therein): 

C 2 (Q) = Co ■ (1 + Ae- Q2/?2 ) • (1 + eQ), (7) 

where Cq is the normalization factor, A is the chaoticity strength factor, 
meaning A = 1 for fully incoherent and A = for fully coherent sources; 
the symbol R is often called as the "correlation radius", and assumed to 
be spherical in this parameterization. The linear term in (JTj) is supposed 
to be account within the long-range correlations outside the region of 
BEC. Note that distribution of bosons can be either far from isotropic, 
usually concentrated in some directions or almost isotropic, and what is 
important that in both cases the particles are under the random chaotic 
interactions caused by other fields in the thermal bath. In the parame- 
terization (j7j) all of these issues are embedded in the random chaoticity 
parameter A. To advocate the formula ([7]) it is assumed: 

a. incoherent average over particle source where A serves to account for: 

- partial coherence, 

- 77 purity; 

b. spherical Gaussian density of particle emission cell (with radius R); 

c. static source which means no time (energy) dependence. 
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However, to enlarge the quantum pattern of particle production pro- 
cess and to avoid the static and undistorted character of particle emitter 
source we have already suggested to use the C^Q) function within QFTp 
accompanying by quantum evolution approach in the form [9]: 

C 2 (Q)^^(iV){l + A 1 (/3)e-^ 2 [l + \ 2 (P)e +Q2R2/2 ]}-F(Q,Ax), (8) 

where £(N) depends on the multiplicity N as £(N) = (N(N - 1))/((N) 2 ). 
The functions \\((3) and A 2 (/3) are the measures of the strength of BEC 
between two photons: Ai(/3) = 7(0;, /3)/(l + a) 2 , and the correction to the 
coherence function in the brackets of Eq. (jHJ) is A2CS) = 2 a/ a/7(w, (3). 
The function j(u, 0) reflects the quantum thermal features of BEC pat- 
tern and is defined as 

^P)=-r{n)= n(a/) , „( w)= n(o;^) = ^- I , u, = 

(9) 

where n(u, (3) is the mean value of quantum numbers for Bose-Einstein 
statistics particles with the energy u in the thermal bath with statisti- 
cal equilibrium at the temperature T = 1/(3. The following condition 
^]-n/(w, (3) = N is evident, where the discrete index / reflects the one- 
particle state /. The function F(Q, Ax) in Eq. flS]) that expresses the 
correlation magnitude as a function of Q and two-particle relative dis- 
tance Ax is a consequence of the Bogolyubov's principle of correlations 
weakening at large distances [12]: F(Q, Ax) = I+77 Q + . . . with normal- 
ization condition F(Q,Ax = 00) = 1; rj is the measure of correlations 
weakening where rf — ^ as A.x - > 00. 

The important parameter a((3) in (jHJ), the measure of chaoticity, sum- 
marizes our knowledge of other than space-time characteristics of the 
particle emitting source, and it varies from to 00 (see [13] for details). 

The mathematical parameter R in formula (jHJ) is the measure of the 
space-time overlap between two photons, and the physical meaning of R 
depends on the fitting of C%(Q) -function. From an experimental point of 
view, three of fitting parameters are regarded as the source size: Rt, Rl 
and i?o, which are transverse, longitudinal and time-like components of 
vector R, respectively. The correlation radius R can be defined through 
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the evaluation of the root-mean-squared momentum Q rms as: 



fo°°d\Q\Q 2 


'C 2 {Q,P)- 


1 


f °°d\Q\ 


~C 2 (Q,P)-1 





C 2 (Q,(3) 



where R and Q rm s(f3) are related to each other by means of (77 - measure 
is neglected) 



R = R{0) 




1 + — J2M 

1 + 4a(/3) V 2 



1/2 



Qrms \P) 



We find the physical range v/372 < (i2 • Qrmsifi)) < V3, where the 
lower bound satisfies the case a — > (no any distortion in the particle 
production domain), while the upper limit is given by the very strong 
influence of chaotic external fields (forces), a — > 00. Actually, the result 
is rather stable in the wide range of variation of a . 

3. It has been emphasized [8,9] that there are two different scale 
parameters in the model considered here. One of them is the so-called 
"correlation radius" R introduced in d7J) and also occurred in (JSJ) . In 
fact, this i?-parameter gives the pure size of the particle emission source 
without the distortion and interaction coming from other fields. The 
another (scale) parameter is the scale L st of the production particle do- 
main where the stochastic, chaotic distortion due to environment (the 
influence of other fields, forces) is enforced. This stochastic scale carries 
the dependence of the particle mass, the a-coherence degree and what is 
very important — the temperature T-dependence. 

One question arises: how can BEC be used to determine the effective 
scale L st and, perhaps, the phase transition? We suppose the changes 
of 77 production region and effects yielding the dynamical variables and 
parameters of BEC due to in-medium distortion. 

Consider the finite system (expanding or squeezing) with the flow of 
two gauge bosons pairs, eg., B^B^- pairs. The Hamiltonian is 



H = - 

° 2 



d x 



2 ^ 2 
+ \VB\ +m 2 B 



which is asymptotically free in the rest frame of undistorted matter with 
the gauge field B^ = (0, B) having the mass m in general case. This 
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Hamiltonian and commutation relations can possess the exact symmetry. 
However, the observed states in real physical environment can not be 
realized in the framework of this symmetry. Hq has to be added by the 
Hamiltonian 

H dis = y — J d 3 xd 3 yB„(x)5F 2 (x-$) B^y), 

which is provided by the distortion due to in- medium effect, in par- 
ticular, because of temperature of the environment. Here, the field 
B^x) propagates in medium with p (momentum) - dependent effective 

(squeezing) frequency Up = yj m 2 + p 2 — 5F^(p) and, consequently, the 
mass m/3 is related to the undistorted (asymptotic) mass m by mp = 
yjm 2 — 5Fp(p). 5F 2 (x) is the non-local formfactor leading to the mod- 
ification of the particle frequency (mass) of undistorted matter. Be- 
cause of the quadratic form of the Hamiltonian Hdi S through the asymp- 
totic Bose - operators of annihilation (creation) a(p)(a + (p)) in B^(x) = 

[ f V, e^.(p) d(p) e~ tpx , one can use the Bogolyubov trans- 

formation 

bj = Uf cif — Vf a^j, — Vja'f — v*j a_/ 

for new operators of annihilation (creation) bf(bi~) in distorted medium. 
In some sense, the last operators correspond to thermalized-distorted 
quasiparticles. The functions u / and v / obey the condition \uf\ 2 — \vf\ 2 = 
1 and can be given in the form Uf — cosh rjf, Vf — sinh rjf, where rjf has to 
be even function of index /. For simplicity, we identify / with momentum 
p and for squeezing frequency one gets r\ v = In u 2 (p)/u^(p). Indeed, the 
Bogolyubov transformation is equivalent to squeezing procedure. 

4. We want to note that in the Higgs rest frame, there is a kinematical 
configuration for two pairs of final lepton momenta (pi,p'j) and (qi,q'j)'- 
Pi ~ p'j and qi ~ q'j. Thus, the final state mimics a two-body final state, 
and if leptons are electrons or even muons, the virtual photons with 
momenta p and p' are nearly on mass-shell, p 2 ~ ~ 4 m 2 , where mi is 
the lepton mass. For the configuration given above, the stochastic scale 
L st has different behavior depending on T. At lower temperatures we 
have 



J st 



16 (2 7r) 3 / 2 e m »/ T 



a(N)M 2 H |l-fe| 2 (m 2 )3/4 T 3/ 2 (l + f ^) 



(10) 
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where the condition n (3 ma > 1 is taken into account for any integer n; 
M H is the mass of the Higgs-boson, 5h = a/1 + (4m M ) 2 /M^; m« stands 
for invariant mass of lepton pair, Am 2 < m\ < (Mh — 2 mi) 2 . 

On the other hand, at higher temperatures when T > nmu, one has 



(4tt) 



3((3)a(N)M 2 I T s \l - 5 H 



|2 



C(3) = ^n- 3 = 1.202. (n; 



n=l 



It turns out that the scale L st defines the range of stochastic forces. This 
effect is given by a(iV)-coherence degree which can be estimated from 
the experiment within the function C 2 (Q) when Q close to zero, C 2 (0), 
at fixed value of mean multiplicity (N) [8,9]: 



a(N) ~ i — . (12) 

C 2 (0)-1 

Note that for C 2 ((5) - function (jHJ) the limit a — > oo yields for fully 
coherent sources while a — > case stands for fully chaotic (incoherent) 
sources. Actually, the increasing of T leads to squeezing of the domain 
of stochastic force influence, and L st (T = T ) = R at some effective tem- 
perature T . The higher temperatures, T > T , satisfy to more squeezing 
effect and at the critical temperature T c the scale L st (T = T c ) takes its 
minimal value. Obviously T c ~ 0(200 GeV) defines the phase transition 
where the chiral symmetry restoration will occur. Since in this phase all 
the masses tend to zero and a — > at T > T c one should expect the sharp 
expansion of the region with L st (T > T c ) — > oo. Our model is consistent 
with the idea of the unification of weak and electromagnetic interactions 
at T > T c , predicted by D.A. Kirzhnitz and A.D. Linde in 1972. In 
addition, at these conditions, there is the analogy with the asymptotic 
free theory: the properties of environment (media) are the same as those 
composed of free particles in the infinite volume (Universe). 

We have also found that increasing of 7 - quantum multiplicity gives 
the raising of a-function that leads to squeezing of L st . 

The qualitative relation between R and L st above mentioned is the 
only one we can emphasize in order to explain the mass and temperature 
dependencies of the source size. 

5. To summarize: one of the main reasons to study BEC (at finite 
T) is the possibility to determine the precision with which the source size 
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parameter and the strength chaoticity parameter can be measured at the 
LHC. It means, for e.g., that to study BEC of two photons one can predict 
the mass of the Higgs-boson. Moreover, such investigations provide an 
opportunity for probing the temperature of the particle production source 
and the details of the external forces chaotic influence. 

In this paper, we propose to calculate the C 2 (Q, (3) - function in which 
the contribution N, T, a, 77 are presented. This differs from the meth- 
ods used in, e.g., LEP experiments based on the approach In fact, 
the latter resembles the traditional way of BEC study, however any qual- 
itative interpretations of R, A, e. The model proposed in this paper is 
expected to be sensitive to the temperature of the environment and the 
external distortion effects. 

We find that the stochastic scale L st decreases with increasing T 
slowly at low temperatures, and it decreases rather abruptly when the 
critical temperature is approached. 

We predict for the first time the spatial size of the source (the corre- 
lation radius R) of two photons originated from the Higgs-boson decays 
in restricted domain at the proper temperature T . The dependence of 
the Higgs-boson mass, the mean multiplicity (N) and the invariant mass 
of the final lepton pair is obtained in the form: 

m/5T 

R ~ (13) 

l/5/\# 2 / 5 M X 1 2 /5 3/10^3/10 v > 

for low values of To < m«, while for higher temperatures one has 

R ~ (14) 

aV5 M J /5 |l-^| 2/5 T 3/5 ' 1 ] 

The correlation radius R decreases as both (N) and mu increase. We 
claim that the experimental measuring of R (in fin) can provide the pre- 
cise estimation of the effective temperature T which is the main thermal 
character in the 7*7* pair emitter source in the proper leptonic decaying 



channel 7*7* — >■ llll with the final lepton energy y kf + mf at given a 
fixed by C 2 (Q = 0) and (N). Actually, T is the true temperature in the 
region of multiparticle production with dimension R = L st , because at 
this temperature it is exactly the creation of two 7 quanta occurred in 
decay of Higgs boson, and these particles obey the criterion of BEC. 

The parameter a can be extracted from the experimental data on the 
two-photon 7*7* BEC using C 2 (Q) with Q being close to zero (see (|T2|) ). 
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This allows one to estimate the Higgs-boson mass, M H , at the temper- 
ature Tq. As a qualitative illustration we present here an estimation of 
Mh assuming To = 100 MeV (the typical temperature of hadronic phase), 
a — 1 (that means the high mean multiplicity case) and mu = 300 MeV. 
The result is rather sensitive to the correlation radius R: Mh = 150 GeV; 
96 GeV and 65 GeV for R = 0.5 fm; 0.6 fm and 0.7 fm, respectively. In 
case of low mean multiplicities (e.g., a = 0.5), we expect the increasing 
of M H as: M H = 213 GeV; 135 GeV and 92 GeV for R = 0.5 fm; 0.6 
fm and 0.7 fm, respectively. 

In the framework of the Longitudinal Center of Mass System (LCMS) 

[14] one can use the average transverse mass (mr) = 0.5(wj9^ + \/Pt 2 ) 



which is applicable to photons with non-zero transverse momentum ^Jp%., 

% — 1,2. In formula ( ...) At is the inverse partial width of the decay 
H — > 77, , and is associated with produced photons through the 
loops (composed of gauge bosons and top-quarks) in the Higgs-boson 
decays, and At represents the strength of electromagnetic interaction 
(see also [15]). 

Because of the fact that (/z + /i~) pairs originated from 7* decay over- 
lap in space and are created in time almost simultaneously, it is natural to 
expect that there are correlation between (u, + u,~) pairs coming from dif- 
ferent 7*'s due to Bose- Einstein interference and the Higgs-boson decays 
in one-loop approximation. These effects may affect the accuracy with 
which the (/i + /i~) pair mass can be measured at the LHC. No systematic 
theoretical treatment of the Bose-Einstein effect within the Higgs decays 
in 7* - pair production has been given so far. In this letter, we presented 
a study of this effect leading to quantum stochastic dependence between 
(fi + fi~) pairs coming from different 7*. 
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Abstract 

It is well understood that the studies of correlations between produced par- 
ticles, the effects of coherence and chaoticity, an estimation of particle emit- 
ting source size play an important role in high energy physics [1]. We mean 
the investigation of the space-time extension or even squeezing of particle 
sources via the multiparticle quantum-statistics correlation. We obtain the 
two-photon correlation function that can provide the space-time information 
about the Higgs-boson source in thermal environment and we argue that such 
an investigation could probe the Higgs-boson mass. 



1. The Large Hadron Collider (LHC) at CERN is at the stage to provide particle 
physicists with treasure of data. These data allowed a precise measurement of many 
important parameters of modern particle physics in order to test their consistency 
and to discover the Higgs-boson. At the same time, little is known about Higgs-gauge 
bosons interplay in particular in theoretical aspects. The effect of Bose-Einstein 
correlations (BEC) is clear and undeniable part of this theory, complicating the 
quantum statistical description of multi-lepton final states. 

Historically, the BEC measurements have concentrated on pion pairs correla- 
tions, however have also been applied to more heavy hadrons, quarks, protons, and 
even gauge bosons, photons, Z-bosons etc. 

It is known that the wave function describing a system of two identical bosons 
should be symmetric under permutation of these bosons. As a consequence, the 
four-momentum differences between two bosons will be smaller than in the world 
where Bose-Einstein statistics would not apply. 

Actually, the direct photon BEC can provide information about the space-time 
distribution of the hot matter prior to freeze-out. However, the BEC with direct 
(primary) photons is faced with difficulties compared to hadron BEC primarily due 
to the small yield of photons emitted directly from the hot region just after particle 



collisions. The main background of photons is produced by decays of final hadrons 
or gauge bosons, or even (pseudo) scalar particles. 

In Higgs searching program at the LHC, we suppose that two photons are pro- 
duced mainly through Higgs-boson decay. From the theoretical point of view the 
decay of Higgs into two photons emerges through the one-loop pattern: the Higgs- 
boson transferred in pair of a quark and antiquark or vector bosons, or even in scalar 
particles. 

In order to reject most of the hadron background, the reconstructed events in 
a spectrometer were required to have a deposit energy (level) of greater than the 
energy well above the minimal energy at which photons emerge from hadrons. Ex- 
perimentally, the 7T° background constitutes a major difficulty which, however, to 
some extent, can be taken care of by measuring the photon pair invariant mass. 

The two-photon correlations in Higgs physics provide a powerful tool to explore 
the Higgs-boson mass estimated via the correlation radius which defines the geo- 
metrical size of the two-photon source. Hence, the two-photon correlation function 
is strongly dependent on the space-time properties of the Higgs decays. 

In the quark-loop scheme, the non-relativistic bremsstrahlung formula for the 
current in quark- antiquark interaction for one photon emission (with four-momentum 
k = (k°, k)) is 

j\k) = ^-p.e x (k)e- k °^, 
m q k {) 

where p is the difference between the spatial momenta of a quark and antiquark with 
the mass m q , £\{k) is the vector of linear polarization of a photon; e means the 
phenomenological parameter which depends on the initial energy. For two sources 
the transition current is 

j\k) = J2e^f n (k). 

71=1 

Index n labels the independent quark-antiquark interactions taking place at dif- 
ferent space-time points y n which are considered to be randomly distributed in 
the space-time volume of the (photon) source. We suppose the gamma-quanta 
are created by some random currents (sources) J u (x) through the Lagrangian den- 
sity Li n t = —gJn{x)A fJ '(x). The neutral weak current of (final) charged leptons is 
J l Ax) = —gl(x)'~f u l(x) with / : e,/j,,r, for which the problem can be solved exactly. 
Both the currents J^x) and Jl(x) exist in a restricted space-time region and they 
are chaotically and randomly disturbed by external fields (forces). The restricted 
domain is characterized by the internal stochastic scale L st , the meaning of which 
is explained in [2,3]. 
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By studying BEC of identical particles, it is possible to determine the time scale 
and spatial region over which particles do not have the interactions. Such a surface is 
called as decoupling one. In fact, for an evolving system such as pp collisions, it is not 
really a surface, since at each time there is a spread out surface due to fluctuations in 
the last interactions, and the shape of this surface evolve even in time. The particle 
source is not approximately constant because of energy-momentum conservation 
constraint. 

Actually, the second order distribution function 



normalized to the product Ni(ki) ■ N 2 (k 2 ) with one-particle distribution functions 
Ni(ki) = (\J Xi (ki)\ 2 ) (i = 1,2) formally defines the probability to find two photons 
with momenta k\ and k 2 issued at yi and y 2 . The crossing momenta has to be taken 
into account. 

In this work, we make an attempt to demonstrate that the problem of properties 
of the genuine interactions can be explored using experimental data which can be 
collected at the LHC. These data can be analyzed through the compared measures 
of some inclusive distributions and final state correlations. 

One of the aims of this paper is to carry out the proposal for the experimental 
measurements of virtual 7*7* pair correlations. 

This exploration is theoretically supported by the quantum field theory model 
approaches [4-9,2,3] at finite temperature, (QFTp), where one of the main parame- 
ters is the temperature of the particle (emitting) source under the random external 
forces (fields) influence. 

We propose that photons do not strongly interact with (produced) medium: 
they carry information about early stage of reaction. Our pragmatic definition is: 
photons are produced not from hadronic decays. Any source of real gamma-quanta 
produces virtual photons with very low mass. If the momentum of the virtual 
gamma-quantum is sufficiently small, the source strength should be small as well. 
The real gamma-quantum can be measured from the virtual yield which is observed 
as low mass of the spectrum for lepton-antilepton pair. 

The main channels are the two-photon production pp — > Higgs — > 7*7* — >• 
2/i~2/i + , 2e~2e + , e~e + /i~/i + , ... in pp collisions. An efficient selection of leptons 
needs to be performed according to the following criteria. First, all leptons were 
required to lie in the pseudorapidity range covered by, e.g., the CMS muon system 
that is, < 2.4. Second, the leptons were required to be unlikely charged in pairs. 

The di-lepton channel is especially promising from the experimental point of 
view, since it is expected that the experimental facilities related for LHC will make 



N 12 (h,k 2 ) = ({^(h) J x *(k 2 )\ 2 ) 
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it possible to record muons of energy in the TeV range with a resolution of about a 
few percent and an efficiency close to 100 %. Moreover, this channel is characterized 
by a maximum signal-to-background ratio in the energy region being considered. 

2. A pair of identical bosons with momenta pi and P2 and the mass m produced 
incoherently from an extended source will have an enhanced probability C 2 (pi, P2) = 
Ni2(pi,P2)/[Ni(pi) ■ N 2 (p2)] to be measured in terms of differential cross section a, 
where 

N 12 {p u p 2 ) = - , n , n 
a alii d\l 2 

to be found close in 4-momentum space 9ft 4 when detected simultaneously, as com- 
pared to if they are detected separately with 



In an experiment, one can account the inclusive density ^2(^1,^2) which describes 
the distribution of two particles in Q (the sub- volume of the phase space) irrespective 
of the presence of any other particles 

, ,11 d 2 n 2 

P2{Pl,P2) 



2! n events dpi dp 2 ' 

where n 2 is the number of particles counted in a phase space domain (pi + dpi,p 2 + 
dp 2 ). The multiplicity N normalizations stand as 



p(p)dp=(N), 



n 



/ P2(pi,P2) dpi dp 2 = (N(N - <5i 2 )), 



where (N) is the averaged number of produced particles. Here, 5 12 = for different 
particles, while 5i 2 = 1 in case of identical ones (coming from the same event). 

On the other hand, the following relation can be used to retrieve the BEC func- 
tion C 2 {Q): 

where N(Q), in general case, is the number of particle pairs (off-shell photons) in 
BEC pattern with 



Q = - V • (pi - P2Y = VM 2 -4m 2 . (2) 
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In definitions ([[]) and (|2J), N re f is the number of particle pairs without BEC and 
p M . = (u>i,Pi) are four-momenta of produced photons (i = 1, 2); M = (pi + j^)^ 
is the invariant mass of the pair of photons. 

An essential problem in two-particle correlations is the estimation of the ref- 
erence distribution N ref (Q) in Eq. CQ). If there are other correlations beside the 
Bose-Einstein effect, the distribution N re f(Q) should be replaced by a reference 
distribution corresponding to the two-particle distribution in a geometry without 
BEC. Hence, the expression ([T]) represents the ratio between the number of 7*7* 
pairs N(Q) in the real world and the reference sample N re f(Q) in the imaginary 
world. Note, that the reference sample can not be directly observed in an experi- 
ment. Different methods are usually applied for the construction of reference samples 
[1] , however all of them have strong restrictions. One of the preferable methods is to 
construct the reference samples directly from data. The 7*7* BEC can be estimated 
for each bin of the photon average transverse momentum p? = \pr x + Pt 2 |/2 as the 
ratio of the distribution of photon pair invariant relative momenta where both pho- 
tons with transverse momenta Pt x and pr 2 were taken from the same event to the 
same distribution but with the photons of the pairs taken from different events. For 
our aim, for the reference sample N re f(Q), it is suitable to use the pairs 7*7* from 
different (mixed) events. 

It is commonly assumed that the maximum of two-particle BEC function C^Q) 
is 2 for pi = f>2 if no any distortion and final state interactions are taking into 
account. 

In general, the shape of C^Q) is model dependent. The most simple form of 
Goldhaber-like parameterization for C^Q) [10,11] is often used for experimental 
data fitting. 

C 2 (Q) = Co ■ (1 + Xe- Q2R2 ) • (1 + eQ), (3) 

where Co is the normalization factor, A is the chaoticity strength factor, meaning 
A = 1 for fully incoherent and A = for fully coherent sources; the symbol R is 
often called as the " correlation radius" , and assumed to be spherical in this param- 
eterization. The linear term in is supposed to be account within the long-range 
correlations outside the region of BEC. Note that the distribution of bosons can be 
either far from isotropic, usually concentrated in some directions, or almost isotropic, 
and what is important that in both cases the particles are under the random chaotic 
interactions caused by other fields in the thermal bath. In the parameterization ([3]) 
all of these issues are embedded in the random chaoticity parameter A. To advocate 
the formula (J3j) it is assumed: 

a. incoherent average over particle source where A serves to account for: 
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- partial coherence, 

- 77 purity; 

b. spherical Gaussian density of particle emission cell (with radius R); 

c. static source which means no time (energy) dependence. 

However, to enlarge the quantum pattern of particle production process and to avoid 
the static and undistorted character of particle emitter source, we have already sug- 
gested to use the C 2 (Q) function within QFTp accompanying by quantum evolution 
approach in the form [2,3]: 

C 2 ( Pl ,p 2 ) ~ £(iV) {1 + Ax(/5) e" A ^ [1 + X 2 (ft) e A ^ 2 ] } , (4) 

where exp(— A p ^) is the smearing smooth dimensionless generalized function with 
Ap» — {.Vi~ V2Y ^fj,u {.Vi~ P2Y ■ ^fiu is the nonlocal structure tensor of the space-time 
size and it defines the domain of emitted photons. £ (N) depends on the multiplicity 
N as f (JV) = (N(N - l))/(iV) 2 . The functions Ai(js) and X 2 (ft) are the measures of 
the strength of BEC between two photons: Xi(ft) = 7(0;, ft)/ (1 + a) 2 , and the correc- 
tion to the coherence function in the brackets of Eq. (j4j) is X 2 (ft) = 2 01/^/7(0;, ft). 
The function j(u},ft) calls for the quantum thermal features of BEC pattern and is 
defined as 

7(c,/3)=7H = nM n{u , y n(u,)=n(u,,ff) = 7 p- T , u, = (5) 

where n(u, ft) is the mean value of quantum numbers for Bose-Einstein statistics 
particles with the energy oj in the thermal bath with statistical equilibrium at the 
temperature T = 1/ft. The following condition rif(co, ft) = N is evident, where 
the discrete index / stands for the one-particle state /. 

The important parameter a(ft) in (jj]), the measure of chaoticity, summarizes our 
knowledge of other than space-time characteristics of the particle emitting source, 
and it varies from to 00 (see [12] for details). 

In terms of time-like Rq, longitudinal Rl and transverse Rt components of the 
space-time size R^, the distribution A p ^ looks like 

A m A pR = (Ap ) 2 Rl + {Ap L f R\ + {Ap T f R 2 T . (6) 

Rq in ([6]) is treated as the measure of the particle emission time, or even it represents 
the interaction strength of outgoing particles. 

Hence, we have introduce a new parameter R^, a four-vector, which defines the 
region of nonvanishing particle density with the space-time extension of the particle 
emission source. Formula (j3J must be understood in the sense that exp(— A p ^) is a 
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distribution that in the limit R — > oo strictly becomes a 5 - function. For practical 
using with ignoring the energy-momentum dependence of a, one has: 



C 2 {Q) ~ £(iV) {l + X x {p) e~ Q2R2 [l + A 2 (/3) e + Q2i?2 / 2 



(7) 



The parameter R in formula flTj) is the measure of the space-time overlap between two 
photons, and the physical meaning of R depends on the fitting of C 2 (Q)-function. 
R can be defined through the evaluation of the root-mean-squared momentum Q rms 
as: 



QrmsiP) = (Q 2 ) 





~C 2 {Q,(3) - 


1 


f °°d\Q\ 


'c 2 {Q,0)-l 





C 2 {Q,P) 



C2(Q,P) 



where R and Q rms (/3) are related to each other by means of 

1/2 



R = R{0) 



1 + 



1 



i 



4a(/3) V 2 



Qrms \P) 



We find the following restricted window < (R ■ Q rm s(/3)) < \/3, where the 

lower bound satisfies to the case a — > (no any distortion in the particle production 
domain), while the upper limit is given by the very strong influence of chaotic 
external fields (forces), a — > oo. The result is rather stable in the wide range of 
variation of a. 

3. It has been emphasized [2,3] that there are two different scale parameters 
in the model considered here. One of them is the so-called "correlation radius" R 
introduced in Q and also presented in (TJJ. In fact, this i?-parameter gives the 
pure size of the particle emission source without the influence of the distortion and 
interaction forces coming from other fields. The other (scale) parameter is the scale 
L st of the production particle domain where the stochastic, chaotic distortion due 
to environment (the influence of other fields, forces) is enforced. This stochastic 
scale carries the dependence of the particle mass, the a-coherence degree and what 
is very important - the temperature T-dependence. 

One question arises: how can BEC be used to determine the effective scale L st 
and, perhaps, the phase transition? We suppose the changes of 7*7* production 
region and effects yielding the dynamical variables and parameters of BEC due to 
in-medium distortion. 
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Consider the finite system (expanding or squeezing) with the flow of two gauge 
bosons pairs, eg., B^B^- pairs. The Hamiltonian is 



(V) 2 + \VB\ 2 + m 2 Bl 



which is asymptotically free in the rest frame of undistorted matter with the field 
Bp = ((p, B) having the mass m in general case. This Hamiltonian and commutation 
relations can possess the exact symmetry. However, the observed states in real 
physical environment can not be realized in the framework of this symmetry. Hq 
has to be added by the Hamiltonian 

H dist = ^ J tfx£yBp{Z)5F}{x - y) B u (y), 

which is provided by the distortion due to in- medium effect, in particular, because 
of temperature of the environment. The field B^(x) propagates in medium with p 

(momentum) - dependent effective (squeezing) frequency oup = ^Jm? + p 2 — 5F 2 (p) 

and, consequently, the mass mp is related to the undistorted (asymptotic) mass m 

by mp = ^m 2 — 5F 2 (p) . Here, 5F 2 (x) is the non-local formfactor leading to the 

modification of the particle frequency (mass) of undistorted matter. Because of the 
quadratic form of the Hamiltonian Hdist through the asymptotic Bose - operators of 
annihilation (creation) a(p)(a + (p)) in 



d 3 p 



one can use the Bogolyubov transformation 

bf = Uf df — Vf alp 6^ = u*f Of — v*f a_/ 

for new operators of annihilation (creation) bf(b~f) in distorted medium. In some 
sense, the last operators correspond to thermalized-distorted quasiparticles. The 

1 2 1 2 

functions Uf and Vf obey the condition \Uf\ — \Vf\ =1 and can be given in the 
form Uf = cosh r]f, Vf = sinh?7j, where f]f has to be even function of index /. 
For simplicity, we identify / with momentum p, and for squeezing frequency one 
gets i] p = In 0J 2 {p)l 'cup 1 (p). Indeed, the Bogolyubov transformation is equivalent to 
squeezing procedure. 

4. In the Higgs-boson rest frame, there is a kinematical configuration for two 
pairs of final lepton momenta (pi,p'j) and (qi,q'j): Pi — p\ and q% ~ q'j. Thus, the 
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final state mimics a two-body final state, and if the leptons being the electrons or 
even muons, the virtual photons with momenta p and p' are nearly on mass-shell, 
p 2 ~ p' ~ 4 m 2 , where m\ is the lepton mass. For the configuration given above, the 
stochastic scale L st has different behavior depending on T. At lower temperatures 
we have 



n 3/2 ^ e 2 mi /T 



A8a(N) m) 1/2 T 3 / 2 



1 + ^, 

16 mi I 



where the condition 2 n (3 mi > 1 is taken into account for any integer n; Mh is the 
mass of the Higgs-boson, mi stands for the lepton mass. 

On the other hand, at higher temperatures, when T > In mi, one has 



TT 2 M 2 H 



_48C(3)a(iV) T 3 mf_ 



C(3) = J>- 3 = 1.202. 



(9) 



n=l 



In case of two real photons correlation one can use the longitudinal stochastic 
scale instead of (jSj) and fl5]), L l °™ 9 {mi = m T /2), with the average transverse mass 

mx = 0.5(y^T + \Jp\ 2 ) in the frame of, e.g., the Longitudinal Center of Mass 

System (LCMS) [13]. 

It turns out that the scale L st defines the range of stochastic forces. This effect is 
given by a(iV)-coherence degree which can be estimated from the experiment within 
the function Ci{Q) as Q close to zero, 6*2(0), at fixed value of mean multiplicity (N): 



1 + Y /2 (n) - C 2 (0) + 7 1/4 (n) A /c* 2 (0)[ 7 1 /2( n )-2]+2 
"(AO = ? — -¥ , (10) 

where 7(71) is defined in (J5J) and C* 2 (0) = C 2 (0)/£(iV). The upper limit of C 2 (0) 
depends on (N) and the quantum thermal features of BEC pattern given by 7(71): 
C*2(0) < £(iV)[l — 7 1 / 2 (n)/2] _1 . This upper limit is restricted by the maximal value 
of 2 in the ideal case as (N) — > 00 and j(n) = 1. 

Note, that for C 2 (Q) - function (j3J), the limit a — > 00 yields for fully coherent 
sources with small (N), while a — > case stands for fully chaotic (incoherent) 
sources as (N) — > 00. Actually, the increasing of T leads to squeezing of the domain 
of stochastic forces influence, and L st (T = T ) = R at some effective temperature 
To. The higher temperatures, T > To, satisfy to more squeezing effect and at the 
critical temperature T c the scale L st (T = T c ) takes its minimal value. Obviously, 
T c ~ 0(200 GeV) defines the phase transition where the chiral symmetry restoration 
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will occur. Since in this phase all the masses tend to zero and a — > at T > T c , one 
should expect the sharp expansion of the region with L st (T > T c ) — > oo. 

Using the relation between Rl and rriT obtained from the Heisenberg uncertainty 
relations [14] Rl(7tit) = C\JhRQ/2f:mT one can estimate Ro within the formula: 

(11) 



R 



c 2 h 



_3C (3) a(N)T*m A T 



where c and h are the speed of light in vacuum and the Planck constant, respectively. 
The expression ffTTj) relates the two-photon emission time Rq with the Higgs-boson 
mass Mh-i the transverse mass of two photons uit and the temperature of the BEC 
domain T > nrnj-- From (TiTJ one obtains Rq ~ 10~ 24 sec and Rq ~ 0.45 • 10~ 24 sec 
for T = 50 GeV and T = 100 GeV, respectively, at M H = 120 GeV and m T = 1 
GeV in the case of two final muon pairs. The confirmation of this should come from 
the measurement of the Higgs-boson lifetime in two-photon decay channel, as well 
as from the BEC data of electroweak interacting particles. 

The qualitative relation between R and L st above mentioned is the only one 
we can emphasize in order to explain the mass and temperature dependencies of 
the source size. The dependence of the stochastic measure of chaoticity a on the 
minimum scale cut L st and Q rm s can be used to define the fit region for different 
Pt- Such a minimum cut on L st introduces a lower cutoff on R. 

There are a number of effects which may give rise to Q ~ correlations and, thus, 
mimic the two-photons BEC coming from Higgs decays. These include i) correlations 
of light hadrons and/or vector bosons misidentified as photons, ii) radiative decays 
of resonances in both pseudoscalar and vector sectors, Hi) collective flow, etc. 

Apparatus or analysis effects which may result in the sensitivity of external 
random forces influence, may be investigated by studying the dependence of the 
correlation functions on the stochastic scale L st . This effect is expected to contribute 
strongly at small L st (or large a and T). 

Our model is consistent with the idea of the unification of weak and electromag- 
netic interactions at T > T c , predicted by Kirzhnitz and Linde in 1972 [15,16]. In 
addition, there is the analogy with the asymptotic free theory: the properties of 
environment (media) are the same as those composed of free particles in the infinite 
volume (Universe). 

5. To summarize: one of the main reasons to study BEC at finite temperature 
is the possibility to determine the precision with which the source size parameter 
and the strength chaoticity parameter(s) can be measured at particle colliders. Such 
investigations provide an opportunity for probing the temperature of the particle 
production source and the details of the external forces chaotic influence. Moreover, 
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one can predict the mass of the Higgs-boson. No systematic theoretical treatment 
of the Bose-Einstein effect within the Higgs decays in 7* - pair production has been 
given so far. 

In this paper, we faced to the model C^Q, 0) - function in which the contribution 
of N, T, a are presented. This differs from the methods used in, e.g., LEP experi- 
ments based on the approach ([3]). In fact, the latter resembles the traditional way of 
BEC study, however any qualitative interpretations of R, A, e have no been clarified. 
The model proposed in this paper is expected to be sensitive to the temperature of 
the environment and to the external distortion effects. 

We find that the stochastic scale L st decreases with increasing T slowly at low 
temperatures, and it decreases rather abruptly when the critical temperature is 
approached. 

We obtain the dependence of the correlation strength functions Ai and A 2 on the 
distance L st and the maximal value of C 2 at minimum invariant momentum Q and 
T. 

We predict for the first time the spatial size of the source (the correlation radius 
R) of two photons originated from the Higgs-boson decay in restricted domain at 
the proper temperature T . The dependence of the Higgs-boson mass, the mean 
multiplicity (N) and the lepton mass is obtained in the form 

R ' 



a i/5 (iV )T 3/10 m; i/1Q 
for low values of T < nmi, while for higher temperatures one has 



R- 



Q,i/5(iv)r 3/5 mf /5 ' 

The correlation radius R increases with heavier Higgs-bosons at large gamma- 
quantum mean multiplicity (N). Actually, the experimental measuring of R (in 
fm) can provide the precise estimation of the effective temperature T which is one 
of the main thermal characters in the 7*7* pair emitter source in the proper leptonic 



decaying channel 7*7* — > Ml with the final lepton energy y kf + mf at given a(N) 
fixed by C 2 (Q = 0) and (N). T is the true temperature in the region of multipar- 
ticle production with dimension R = L st , because at this temperature it is exactly 
the creation of two 7 quanta occurred in decay of Higgs-boson, and these particles 
obey the criterion of BEC. 

The parameter a can be extracted from the experimental data on the two-photon 
7*7* BEC using C 2 (Q) with Q being close to zero (see ( ITU]) ). This allows one 
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to estimate the Higgs-boson mass, M H , at the temperature T Q . As a qualitative 
illustration we present here an estimation of M H assuming T = 50, 100 GeV, a = 10 
% and mi = m M = 105 MeV. The result is rather sensitive to the correlation radius 
R and T : M H = 167 GeV at T = 50 GeV, M H = 477 GeV at T = 100 GeV for 
R = 1.0 fm; M H = 68 GeV at T = 50 GeV, M H = 194 GeV at T = 100 GeV for 
R = 0.7 fm- 

Because of the fact that (fi + ^) pairs originated from 7* decay overlap in space 
and are created in time almost simultaneously, it is natural to expect that there are 
correlation between pairs coming from different 7*'s due to Bose-Einstein 

interference. These effects may also affect the accuracy with which the pair 
mass can be measured at the LHC. 
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